Controlling the efficiency of spin injection into graphene by carrier drift 
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Electrical spin injection from ferromagnetic metals into graphene is hindered by the impedance 
mismatch between the two materials. This problem can be reduced by the introduction of a thin 
tunnel barrier at the interface. We present room temperature non-local spin valve measurements 
in cobalt/aluminum-oxide/graphene structures with an injection efficiency as high as 25%, where 
electrical contact is achieved through relatively transparent pinholes in the oxide. This value is 
further enhanced to 43% by applying a DC current bias on the injector electrodes, that causes carrier 
drift away from the contact. A reverse bias reduces the AC spin valve signal to zero or negative 
values. We introduce a model that quantitatively predicts the behavior of the spin accumulation in 
the graphene under such circumstances, showing a good agreement with our measurements. 

PACS numbers: 72.25.Hg, 73.63.-b 
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The predictions of a long spin relaxation time in 
graphene [l| , together with the availability of micrometer 
sized graphene flakes 0] fueled a number of experimen- 
tal studies on graphene spin valve type devices. Since 
2006, several successful spin injection experiments were 
reported in field effect transistor geometries of lateral 
dimensions of a few micrometers H, 0, [E S S B @j 
all relying on electrical spin injection and detection us- 
ing ferromagnetic metal electrodes. In earlier experi- 
ments Q we have determined the spin polarization of our 
ferromagnet/aluminum-oxide/graphene contacts as 10% 
(room temperature experiments), strongly depending on 
the contact resistances i.e. the properties of the tunnel 
barrier. We have also presented a method to manipulate 
the spin transport through graphene from injector to de- 
tector by means of carrier drift, under the influence of a 
DC electric field 

At this point the question arises, what determines the 
efficiency of spin injection from the ferromagnet into 
graphene. In the limit of clean metal on graphene, exper- 
iments show that an ohmic contact is formed at the inter- 
face, and spin injection is determined by the spin selec- 
tive resistivity of the ferromagnet. Due to the impedance 
mismatch between the two materials [111, [13 ) this leads 
to a very inefficient injection of the spin polarized cur- 
rent into the graphene. In case of a tunnel barrier, on the 
other hand, the spin dependent tunneling takes the role 
of the spin-dependent contact resistance, and impedance 
mismatch can be reduced jl3|]. This way we can define 
an intrinsic spin polarization of a contact P that is deter- 
mined by the nature of the interface/tunnel barrier and 
the ferromagnet, and an effective injection (or detection) 
efficiency Pmj/det that takes into account the presence of 
graphene and is lower due to the impedance mismatch. 

A study of the interplay of mechanisms behind electri- 
cal spin detection in lateral Fe/GaAs structures (Schot- 
tky tunnel barriers) was done recently by Crooker et 
al.[l3l- The authors found that the sensitivity of such a 
spin detection scheme can be tuned by an electrical bias 
on the interface. The effect was explained by the bias 



dependence of the tunneling spin polarization (interface 
effect) as well as the bias dependence of the spin trans- 
port in the GaAs before the spins reach the detection 
point (bulk effect). 

Here we present the manipulation of the effective 
spin injection efficiency at electrically biased ferromag- 
net /tunnel barrier/graphene interfaces, where the tunnel 
barrier is below 1 nm thick with (probably) relatively 
transparent pinholes. The four-terminal spin valve de- 
vices built for this purpose, illustrated in Fig. [Tla), are 
similar to the ones investigated in Ref. 3 and see 
these references for a detailed description of the fabrica- 
tion procedure. In the non-local measurement geometry 
we employ the current injection circuit (F3, F4, F5) is 
separated from the voltage probes (Fl, F2) decoupling 
the charge current from the spin current, a technique 
that allows detection of a voltage difference that can only 
be attributed to a spin diffusion through the graphene 
laycr[15j . The spin injector/detector electrodes are 50 nm 
thick Co strips of widths (left to right) 800, 250, 90, 140 
and 350 nm. The different magnetic shape anisotropics 
yield switching fields from below 20 mT (Fl, widest) up 
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FIG. 1: (color online) (a) The sample layout and non-local 
measurement geometry, injection schemes A (black) and B 
(red). The direction of arrows indicate positive bias, (b) 
Illustration of the carrier drift on a length scale L from a pin- 
hole of width w. Under DC bias, a strong local electric field 
induces drift of carriers in the graphene (see arrows) facili- 
tating or blocking the spin injection from the ferromagnetic 
electrode F3 through the pinhole. 
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to 85 mT (F3, narrowest). The Co contacts patterned 
by electron beam lithography are separated from the 
graphene by a 0.8 nm thick AI2O3 layer; the distances 
between them are shown on the figure. The contact re- 
sistances probed by three-point electrical measurements 
were between 50 kft and 90 kil. The fact that the con- 
tact resistances do not scale with the contact areas, as 
well as the granular morphology of the sub-nanometer 
aluminum-oxide layer indicates an electrical contact me- 
diated by a number of pinholes. The spin injection and 
detection method was based on a standard low-frequency 
AC lock-in technique with a current of 0.5 /iA RMS. In 
addition to the AC injection current, a DC current bias (- 
5 to 5 ^A) was applied on the current injector electrodes 
F3, F4 and F5, in a geometry explained below. 

In the first set of experiments we used contacts F3 and 
F4 as AC-I-DC current injectors, and contacts F2 and Fl 
as the non-local voltage probe (scheme A on Fig. [Ija)). 
As a consistency check, we performed a second set of 
experiments on the sample, where everything was kept 
the same except the current injectors, that we shifted 
to the electrodes F4 and F5, leaving electrode F3 float- 
ing (scheme B). This resulted in an opposite DC bias on 
the electrode F4 in comparison with scheme A. Further- 
more, a second sample was manufactured with the same 
spin valve geometry but contact resistances in the 1 - 
10 kS7 range, where we have performed similar measure- 
ments. The results support the data we present in this 
manuscript, indicating good reproducibility of the effect. 

In these experiments we swept the magnetic field 
aligned parallel to the electrodes from -150 to -fl50 
mT, while monitoring the non-local resistance defined as 
= Vac /I AC- This was repeated for a number of se- 
lected bias currents Idc- The two measurement sets for 
schemes A and B are shown on Fig. [2l Every step in the 
spin valve measurements, as indicated for the case of the 
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FIG. 2: (color online) Non-local spin valve measurements at 
DC current biases -5, -2.5 , 0, 2.5 and 5 ^A [panels (i) to (v)], 
for contact scheme A (left column) and B (right column). In 
scheme A, the magnetic switching of contacts F2, F3 and F4 
is visible (also shown by black arrows). 



zero bias curves (iii), corresponds to the magnetization 
switch of an electrode. Note, that electrodes Fl and F5 
were at a distance approximately 3 times greater than 
the spin diffusion length in graphene (As/ — 2/^m 0]), 
therefore their contribution to the measured resistance 
was too weak to appear in the plotted data. By electro- 
static measurements the graphene was determined to be 
p-type, carrier density being in the 10^^ range (at 
zero apphed gate voltage). 

Examining the zero-bias measurement in scheme A 
(panel iii left), we can follow the evolution of the re- 
sistance while sweeping the magnetic field from zero to 
-1-150 mT. We start with all four electrodes Fl-4 aligned 
parallel ("up" direction) in their magnetization. In this 
scheme, electrode F3 injects spin-up polarized current 
and creates a spin imbalance in the graphene. Electrode 
F4 extracts spin-up, i.e. creates an opposite spin imbal- 
ance. The spin imbalance diffuses through the graphene 
and arrives to the detector F2 (sensitive to the spin-up 
channel). Since F3 is much closer to the detector than 
F4, its effect on the detector will be stronger and we 
measure a positive non-local resistance of approximately 
2517. When the magnetic field reaches the value of 20 
mT, we see a step in the resistance that we associate 
with the switching of detector F2 to the "down" direc- 
tion, being now antiparallel with the injector F3. F2 
probes thus the " down" spin channel, yielding a negative 
resistance level. At approximately 50 mT, the injector 
electrode F4 switches and another step appears in the 
resistance. The two injectors F3 and F4 are now an- 
tiparallel oriented, which means they both inject spin up 
carriers. The switching of F4 therefore further increases 
the signal on the detector electrodes, this becoming more 
negative. Finally, at 85 mT field the injector F3 switches 
to the "down" state. All four electrodes being parallel 
again, the measured resistance shows the initial value of 
+25f7. 

On the right panel (iii) of Fig. [21 in the injection scheme 
B, we see only two steps, at magnetic fields of 20 and 50 
mT. Electrode F3 being not connected, it does not con- 
tribute to the signal. The four vertical arrows represent 
the magnetic orientation of Fl, F2, F4 and F5. The 
steps associated with the switching of electrodes F2 and 
F4 happen at the same magnetic fields both in scheme A 
and B, indicating a consistent behavior of the spin valve. 

The non-local resistance levels we measure depend on 
the relative magnetic orientation of the electrodes con- 
trolled by the magnetic field, on the distance between 
injectors and detectors fixed by design, and on their spin 
injection/detection efficiencies. The efficiencies are cor- 
related to the spin polarization of the contacts between 
graphene and injectors Bmj and, respectively graphene 
and detectors Pdet- These can be calculated from the 
spin valve measurement using the relation from Ref. 

^nl = ^^1^^^-P(-^-/^^/)' 
where the non-local resistance R^^ is given by the mea- 
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FIG. 3: (color online) (a)The resistance steps associated with 
the switching of electrodes F4,F3 in the schemes A and B, 
plotted as a function of the applied DC current bias; (b) sim- 
ulated with the drift model for two different drift lengths L. 



surements, \sf — 2fim is the spin diffusion length in 
graphene 3|, Wg = 500nm is the width of the graphene 
channel, a — ~ 1.1 • 10~^r2~^ is its measured con- 
ductivity and Lsv is the injector-detector distance (F3 to 
F2 and F4 to F2, respectively). For our case, the spin 
polarization this formula yields is 25%. 

Applying an electric bias on the injector electrodes F3 
and F4 (scheme A) and respectively, F4 and F5 (scheme 
B) changes the picture. On Fig. [2] panels (i), (ii), (iv) 
and (v) we plotted the spin valve measurements taken at 
DC bias currents of -5, -2.5, -1-2.5 and -1-5 /xA. Through 
the measurement series (i) - (v) in scheme A we can fol- 
low the behavior of the resistance steps (thus the spin 
injection/detection efficiencies) associated with the elec- 
trodes F2, F3 and F4 under the applied biases. In case of 
injector electrode F3, a positive bias (panels iv, v) does 
not seem to change much in the stepsize; it stays -1-75 fl. 
However, when the bias is reversed, the resistance step 
is gradually reduced to zero. This indicates, that under 
negative bias, the spin injection of this contact becomes 
very inefficient. The behavior of the injector electrode F4 
is consistent with its opposite bias compared to F3. The 
resistance step associated with it is increasing from -100 
O to -1-10 when we change the bias from -5 /lA to +5 
/iA. Note, that the steps associated with the efficiency of 
F4 are generally lower, due to the larger injector-detector 
distance, but they follow approximately the same trend 
as in the case of F3. The sensitivity of detector F2 is 
unchanged during the measurements (no bias applied on 
the detectors); the resistance steps associated with its 
switching are simply equal to the difference of the two 
steps caused by F3 and F4. For a better overview, on 
Fig. [3] (a) we have summarized these observations, plot- 
ting the behavior of the resistance steps against the cur- 
rent bias. The measurements done in scheme B deliver 
an additional curve for F4; this electrode is now under 
positive bias, and the sign of both the spin valve signal 
as well as Idc is reversed. 

These experiments indicate, that biasing the injector 
electrodes yields a dramatic change in their spin injec- 
tion efficiencies, enhancing the spin valve signal to a 
saturation value or suppressing it completely. Applying 
Equation (1) for the spin valve signal we measured in 
scheme B, at -1-5 fiA bias we calculate a spin polarization 



Pinj = 43% for the contact between F4 and graphene. 
To do this, we keep the spin polarization of the unbiased 
detector/graphene contact Pdet at the original 25%. In 
case of maximum reverse bias on the other hand, the AC 
measurements show a reversed spin valve behavior (see 
panels (i) on Fig.[2|) that will be addressed later. 

The physics behind spin injection through an inhomo- 
geneous tunnel barrier under the action of a DC elec- 
tric field can be modeled and understood considering a 
strong local drift of charge carriers. This happens in the 
graphene region directly around the injection points (the 
pinholes) on a characteristic length scale L. The width of 
the pinholes w in the AI2O3 barrier is another important 
parameter that defines the electric field E possible to be 
generated by a current I^c sent through the contact. An 
illustration of the idea is included in Fig. [1] (b). 

A qualitative picture of the effect can be created in 
the following way. From injection to detection, there are 
three electrically coupled regions of different spin trans- 
port: injection from the ferromagnet through the tunnel 
barrier via pinholes, spin transport through the graphene 
in the immediate vicinity of the injection point (drift re- 
gion) and finally, diffusion/relaxation towards the detec- 
tion point. In the drift region the DC bias gives rise to the 
electric field E. This yields a drift-diffusion type of trans- 
port as described in Ref. [l3] , however in this case on the 
short length scale L. The carriers in the graphene drifting 
away from the injection point (in case of positive DC bias 
and p-type graphene) reduce the backflow of spins and 
thus facilitate further injection of spin polarized current 
through the pinhole(s). The upper limit of the effect is 
a measurement of the intrinsic spin polarization P possi- 
ble to inject from the ferromagnetic electrodes (large E), 
when impedance mismatch is eliminated. On the other 
hand, an opposite electric field polarity (or carrier type) 
will result in a carrier drift towards the injection point 
and therefore, it will keep the local spin polarization 
high, reducing the efficiency of the spin-polarized injec- 
tion. An increased negative bias enhances the impedance 
mismatch to a point when the drift starts to dominate 
the spin transport. In our AC measurements, this yields 
a negative differential resistance (lower spin injection for 
stronger electric fields) and thus the reversed AC spin 
valve signals on Fig. [51 Finally, above a threshold elec- 
tric field value the strong drift effect prevents detection 
of any AC spin signal at the electrodes F2,F1. 

To quantitatively model the transport in the drift re- 
gion, we consider a one-dimensional drift-diffusion of the 
spins along the x-axis. We write the spin accumulation 
Us along the axis x {x = represents the injection point) 
as described in Refs. [l^, [ill in the form of 



n,.(x) = Aexp(+-^) -t-Bexp(--^), (2) 

A-L A — 



where X± are the up/downstream spin transport length 
scales as in fl^ . The spin current density flowing through 
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the region due to diffusion and drift can be written as 

3,{x)^-D'^:^^^^+VD-n,{x). (3) 
ax 

where vd represents the drift velocity. 

The coupHng at the edges of the drift region yields 
two boundary conditions. Injection from the Co elec- 
trode is modeled by a spin-polarized current source 
Is —PI- where P represents the intrinsic spin po- 
larization without any impedance mismatch - parallel 
with a spin flip resistor Rg of dimension ■ s. This 
contains the contact resistance Rc — 50 kfi and the 
impedance mismatch between the graphene and the fer- 
romagnet: Rg = ARcWgN2De'^ , where Wg = 500nm is 
the graphene width and N2D is the 2-dimensional den- 
sity of states in graphene. The relation iV2De^ — D/Rg 
allows us to determine Rs from the measured values of 
the graphene sheet resistance Rg ~ 875f2 and diffusion 
constant D ^ 0.02m2/s. The spin current density enter- 
ing the drift region is therefore the spin current of the 
source Is minus the relaxation through the contact rep- 
resented by Rs- 

The spin current density exiting the drift region is cou- 
pled to the other end of the device, where spin transport 
is governed by diffusion and relaxation, and we model 
it with another spin flip resistor Rout = ^sf I'^D. Here 
the factor 2 is introduced since spin relaxation can take 
place on both sides of the contact. These considerations 
allow us to calculate the spin accumulation at the detec- 
tor electrode, as a function of the electric field E and the 
drift length L. Knowing the spin polarization of the de- 
tector electrode P, we obtain the AC voltage difference 
Vac on it that is due to the spin transport. On Fig.[3](b) 
we plotted this voltage difference in form of a non-local 
resistance Vac /Iac against the electric field for two dif- 
ferent drift lengths, using the experimental parameters of 
the injection scheme B. A linear correspondence between 
the X-axis of Fig. EJ^a, experiment) and (b, simulation) is 
given by the formula E — I^c ■ Rl I L = Idc • Rg/w, where 
Rl is the total resistance of the graphene along the drift 



region L and w is the pinhole size. In this case, we refer 
to the collective effect of all pinholes in the tunnel bar- 
rier that contribute to the injection process. Comparing 
the simulation to the measurements and considering the 
length of the modeled drift region L to be around 100 
nm, we note that the electric fields necessary to obtain 
good accordance are in the ±10^ V/ni range. This results 
in a drift velocity ■yu = M ' — 0-25 • 10^ m/s that equals 
approximately the quarter of the Fermi velocity. This is 
only possible to achieve with w; < 5 nm. Considering, 
that we probably have more than one pinhole, the indi- 
vidual average size of them is at the nanometer range or 
below. This is consistent with the fact, that spin trans- 
port under the Co contacts is possible (the graphene is 
not shorted by large contact areas). 

In the presented measurements, the entire graphene 
channel was p-type. We have done similar sets of mea- 
surements where the graphene between the electrodes 
was n-type and, respectively in the vicinity of the Dirac 
neutrality point, using electrostatic gating. Aside of a 
scaling of the signal with the carrier density, the results 
were consistently the same in all three regimes. This 
means that gating has no influence on the spin injection 
process, only on the diffusive transport in the graphene 
channel between the contacts [l7j . 

In conclusion, we demonstrated electrical spin injec- 
tion into graphene at room temperature with a high effi- 
ciency of 25%, controllable by electrical bias on the con- 
tacts. Non-local signals up to lOOfi were obtained this 
way. The results were explained by injection through 
nanometer sized pinholes in the tunnel barrier between 
the metal and graphene followed by strong local carrier 
drift. We do not fully understand yet the nature of our 
AI2O3 barrier. The contact resistances we measured in- 
dicate the presence of relatively transparent (though not 
metallic) pinholes; the morphology of the oxide layer has 
yet to be verified by surface characterization experiments. 

This work was financed by the Zernike Institute for 
Advanced Materials, NanoNed, NWO and FOM. 
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